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It has been proposed that the radio-frequency electric-dipole (El) transition between two nearly 
degenerate opposite-parity states in atomic dysprosium should be highly sensitive to possible tem- 
poral variation of the fine structure constant (a) [V. A. Dzuba, V. V. Flambaum, and J. K. Webb, 
Phys. Rev. A 59, 230 (1999)]. We analyze here an experimental realization of the proposed search 
in progress in our laboratory, which involves monitoring the El transition frequency over a period 
of time using direct frequency counting techniques. We estimate that a statistical sensitivity of 
\a/a\ ~ 10~ 18 /yr may be achieved and discuss possible systematic effects that may limit such a 
measurement. 

PACS numbers: PACS. 06.20.Jr, 32.30.Bv 



I. INTRODUCTION 

Variation of the fundamental constants of nature would 
signify new physics beyond the Standard Model, as dis- 
cussed in a recent review 0. Various theories con- 
structed to unify gravity with the other forces allow or 
necessitate such a variation 0EJ[3 Of recent interest 
is the astrophysical evidence for a variation of the fine 
structure constant a. From an analysis of quasar ab- 
sorption spectra Q over the redshift range .5 < z < 3.5, 
a 4cr deviation of Aa/a = (-0.72 ± 0.18) x 10~ 5 from 
zero was reported. Although these data hint at nonlin- 
ear dependence of a with time, for simplicity we assume 
a linear shift over ~ 10 10 years, which corresponds to a 
temporal variation of a/a = (7.2 ±1.8) x 10~ 16 /yr. The 
current best terrestrial limit, over a much shorter time 
scale of 2 billion years, is \a/a\ < 10~ 18 /yr 0, H ||], 
which comes from an analysis of geophysical data ob- 
tained from a natural fission reactor at Oklo (Gabon) 
which operated 1.8 x 10 9 years ago. Observational mea- 
surements like these, however, are subject to numerous 
assumptions which tend to complicate the interpretation. 
For example, questions have been recently raised regard- 
ing the reliability of the Oklo analysis |lCj . 

Laboratory searches have numerous advantages, but 
have, thus far, placed weaker limits on a. For example, 
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a comparison between H-maser and Hg + gave a/ a < 
3.7xl0~ 14 /yr 11] while a similar comparison between Rb 
and Cs microwave clocks yielded \a/a\ < 1.6 x 10~ 15 /yr 
[T^ . A limit of | d/a | < 1.2 x 10~ 15 /yr was obtained from 
a com par ison of a Hg + optical clock to a Cs microwave 
clock |13| using a frequency comb. 

It has been suggested ^| that the electric-dipole 
(El) transition between two nearly degenerate opposite- 
parity states in atomic dysprosium (Dy; Z=66) should be 
highly sensitive to variations in a. Indeed, a recent cal- 
culation 

M supports this conclusion. An experimental 
search utilizing these states is currently underway and is 
discussed here. We provide an analysis of possible sys- 
tematic effects and show that this search could ultimately 
reach a sensitivity of \a/a\ ~ 10~ 18 /yr. 

II. VARIATION OF ALPHA IN DYSPROSIUM 

Tests of variation of a in atomic systems rely upon the 
fact that rclativistic corrections depend differently on a 
for different energy levels. The total energy of a level can 
be written as 

E = E + q(a 2 /a 2 - 1) + 0(a 3 ), (1) 

where Eq is the present-day energy, ao is the present-day 
value of the fine structure constant, and q is a coefficient 
which determines the sensitivity to variations of a. This 
coefficient mainly depends upon the electronic configura- 
tion of the level. A recent calculation 16], utilizing rela- 
tivistic Hartree-Fock and configuration interaction meth- 
ods, found values of q for two nearly degenerate opposite- 
parity states in atomic dysprosium that are both large 
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FIG. 1: Partial energy-level diagram of dysprosium showing the present scheme to populate level B and detect the population 
of level A. Solid lines: absorption. Dashed lines: spontaneous emission. 



and of opposite sign. For the even-parity state (desig- 
nated as A), qA/hc — 6008 cm -1 while, for the odd- 
parity state (designated as B), qs/hc = —23708 cm -1 . 
The time variation of the transition frequency between 
levels A and B is (for a w arj) 



v = 2 



qB - qA 



a/a 2 x 10 15 Hz a/a. 



(2) 



In other words, \a/a\ = 10~ 15 /yr, implies \i>\ — 2 Hz/yr. 

The statistical uncertainty of determining the central 
frequency from a resonance lineshape is 5v ~ j/\/~N, 
where N is the number of counts and 7 is the tran- 
sition width. In the dysprosium system, 7 = 20 kHz 
is determined mainly by the lifetime of state A, ta — 
7.9 (is (jb > 200 [is) 01- A reasonable counting 
rate of 10 9 s -1 , corresponds to a statistical sensitivity 
to i>o of ~ 0.6 t -1 / 2 Hz^fs where r is the integration 
time in seconds. After an integration time of one day, 
8u ^ 2 mHz, thus allowing for a statistical sensitivity of 
\a/a\ ~ 10 _18 /yr for two measurements separated by a 
year's time. 

It should be noted that the sensitivity to a variation 
in a is not proportional to 5v/v. In fact, if v is suffi- 
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ciently small, the uncertainty in its determination is no 
longer limited by the relative uncertainty of the refer- 
ence clock frequency {y clock) ^ but rather by other exper- 
imental uncertainties. For example, a modest Cs clock 
provides an absolute accuracy of 10 -12 . In the Dy tran- 
sitions considered here, typical frequencies are ~ 1 GHz 
or smaller. Thus the clock uncertainty will only become 
an issue when other uncertainties can be reduced to the 
level of 
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Hz, 



(3) 



and even then can be overcome in a straightforward way 
by using a better reference clock. 

The relatively modest requirement for the reference 
clock also means that the Dy experiment is insensitive 
to a variation of the clock's frequency due to a change 
in fundamental constants. For example, if a varies by 
10~ 15 , the fractional change in the Cs clock frequency 
will be 2 x 10~ 15 , as the hyperfine transition frequency 
employed in the clock is proportional to a 2 . This vari- 
ation is negligible in comparison with the clock stability 
oflO" 12 . 
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III. EXPERIMENTAL TECHNIQUE 

A. Overview 

The experimental search for a variation of the El tran- 
sition frequency between the two opposite-parity states 



will proceed as follows. As shown in Fig. atoms are 
populated |lg to the longer-lived odd-parity state B 
(jb > 200 fj,s 01) v i a three transitions. The first two 
transitions require light at 833 nm and 669 nm, respec- 
tively. The third transition involves spontaneous decay 
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FIG. 2: Hyperfine structure of of A and B levels of 163 Dy. 
Zero energy is chosen for the lowest hyperfine component. 



with a branching ratio of 30% and the emission of 1397- 
nm light. Once in state B, atoms are transferred to state 
A with an rf electric field, whose frequency is referenced 
to a commercial Cs frequency standard. A resonance 
lineshape is attained by scanning the rf frequency and 
monitoring 564-nm light from the second step of fluores- 
cence from state A. 



B. rf transitions 

Because the energy separation between the nearly de- 
generate levels is on the order of hyperfine splittings and 
isotope shift energies (Fig. fJJ), and because dysprosium 
has seven stable isotopes, there are many choices of rf 
frequencies. Table[I]shows rf transitions with frequencies 
below 2 GHz calculated from measured hyperfine con- 
stants and isotope shifts 01 • The smallest transition 
frequency (3.1 MHz) occurs for the F — 10.5 components 
of 163 Dy (the same transition is used in a search for parity 
nonconservation (PNC) 20] ). A low- frequency transition 
that offers a higher counting rate is the 235 MHz transi- 
tion in 162 Dy. This is due to a large isotopic abundance 
and the fact that there is no hyperfine splitting to dilute 
the atomic population. As discussed below, the choice of 
rf transition is also important in regard to sensitivity to 
systematic effects. 



C. Apparatus 

We have studied states A and B extensively as a sys- 
tem to measure atomic PNC effects |20|. Although our 
current apparatus is optimized for a PNC experiment, it 
is suitable for a measurement of a with only minor mod- 
ifications. We describe this system here in order to make 
a realistic evaluation of a possible experiment. 

An atomic beam is produced by an effusive oven op- 
erating at 1500 K . The atoms pass through collimators 



TABLE I: Calculated El transition frequencies using the hy- 
perfine constants and isotope shifts from Ref. |l7j|. ||dF|| and 
||dj|| = 1.5(1) x 10~ 2 ea « 19 kHz/(V/cm) are reduced ma- 
trix elements. The last three columns list the isotope mass 
numbers (and abundances) and total angular momenta of the 
states involved in the transition. 
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and approach the interaction region (Fig. |2Jl where the 
desired electric and magnetic fields are produced. In- 
side the interaction region, the atoms encounter the laser 
beams used in the population scheme. The rf electric 
field is formed between two wire grids which are used in 
order to minimize surface area and thus stray charge ac- 
cumulation. Typical applied electric-field amplitudes are 
~ 5 V/cm. If necessary, a magnetic field of up to ~ 4 G 
can be applied parallel to the electric field. It is produced 
by eight turns of gold-plated copper wire surrounding the 
interaction region. The entire interaction region is placed 
inside a magnetic shield. As atoms decay from state A, 
564-nm light from the second step of fluorescence (Fig.^) 
is collected by a light pipe and detected by a photomul- 
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TABLE II: Estimated sizes of systematic effects thus far con- 
sidered. 
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FIG. 3: Side view of interaction region currently optimized 
for PNC search. 
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FIG. 4: Zeeman-level-crossing resonances for 163 Dy. State B 
is populated as shown in Fig. Q in the presence of an electric 
field of 0.4 V/cm. The population of state A, resulting from 
electric field-induced mixing of states A and B, is monitored 
by observing 564-nm fluorescence (Fig. 0. The solid line 
represents the best curve fit using the appropriate initial level 
separation [rH ] and g- values |Tflj |. 



tiplicr tube. The extent of the light pipe area is shown 
by the dashed box in Fig. |3J 

An example of a resonance lineshape attained with this 
apparatus is shown in Fig. 0] Here, a magnetic field is 
scanned in the presence of a dc electric field revealing 
Zeeman-level-crossing resonances for 163 Dy. For the mea- 
surement of a, the signal-to-noise will be much greater 
as individual Zeeman resonances will not be resolved. 

The interaction region can be improved by making it 
shorter, as most atoms decay within twice the lifetime 
of state A (ta = 7.9 /is). Given a mean atomic velocity 
of 5 x 10 cm/s, an appropriate length will be ~ 8 mm. 
Furthermore, a shorter interaction region will allow for 
improved light collection efficiency, better suppression of 
background oven light, and better control over E- and 
B-fields. 
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D. Lineshape 

The natural lineshape is a Lorentzian with a width 
(7 ss 20 kHz) determined mainly by the lifetime of state 
A (ta = 7.9 /is). To a first approximation, both transit- 
time and Doppler broadening change the lineshape in a 
symmetrical fashion. For an interaction length ~ 10 cm 
and a mean atomic velocity of v — 5 x 10 4 cm/s, the 
transit-time width is ~ 5 kHz. The Doppler shift (for an 
atom moving along the propagation direction of a free 
electromagnetic wave of frequency u) goes as ~ is(v/c), 
which for a 1 GHz transition is ~ 2 kHz. Thus, both 
transit-time and Doppler broadening contributions to 
the width are much smaller than the natural linewidth. 
Moreover, note that symmetrical broadening of the reso- 
nance does not change the central frequency. In addition, 
although an asymmetrical resonance can cause an appar- 
ent shift of the central frequency, this is not important 
for the search for variation of a as long as the asymmetry 
does not change between measurements. 

In the following section, we discuss mechanisms lead- 
ing to possible lineshape asymmetry and estimate the 
corresponding shifts. 



IV. SYSTEMATIC EFFECTS 



High statistical sensitivity implies that our technique 
will likely be limited by how well we can control system- 
atic effects. Below we analyze possible systematics and 
give estimates of how well various parameters of the ex- 
periment need to be controlled in order to achieve a sen- 
sitivity of \a/a\ ~ 10 _18 /yr. A summary of systematic 
effects and their estimated sizes is given in Table UTI 
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A. dc Stark and quadrupole shifts 

Variations of a stray dc electric field (Ed c ) can give 
rise to time-varying frequency shifts for the transition 
frequency between states m and n: 



5v dc = 



rft P 2 

a mj- Cj dc 



A 



(4) 



where d m j and A m j = E„ 



- Ej are the dipole matrix ele- 
ment and energy separation, respectively, between states 
m and j. The sum in Eq. (0J) is taken over all states, 
including hyperfme levels in the case of an odd isotope. 

The reduced dipole matrix element for the B — > A tran- 
sition is \\dj\\ = 1.5(1) xlO" 2 ea « 19 kHz/(V/cm) [l7|. 
For large J, the maximum z projection of the dipole mo- 
ment w ||dj||/\/2J = 4kHz/(V/cm). Matrix elements 
connecting all levels with the same configurations as A 
and B have similar values pl|. Thus, between levels 
within 2000 cm -1 of the A and B levels, dipole matrix el- 
ements are relatively small, which reduces the sensitivity 
to systematics related to stray electric fields. 

In an earlier PNC search, we reported stray dc E-fields, 
presumably due to stray charge accumulation on the elec- 
trode surfaces, that varied on time-scales from hours to 
months, and had a typical magnitude of 50 mV/cm pfj) . 
Using this value for the field and d ~ 4 kHz/ (V/cm), we 
estimate the shift for transition frequencies in the range 
3-1000 MHz to be - 10~ 4 - 10~ 2 Hz. The stray electric 
field can also be measured at a few mV/cm level using the 
atoms themselves |20j and cancelled by the application 
of an external electric field. 

An important systematic in Hg + optical clock experi- 
ments |22| is the electric quadrupole shift due to a stray- 
field gradient (VE). Based upon the size of the inter- 
action region and the homogeneity of the electric field 
for the Dy beam experiment (10~ 3 ) 20], a conservative 
estimate gives VE ~ 10 mV/cm 2 . Assuming a typical 
a quadrupole moment Q ~ 1 ea 2 , leads to a quadrupole 
shift estimated to be < 10~ 5 Hz which is negligible. 



B. ac Stark shift 

Fluctuations of the rf electric-field amplitude (E) can 
also lead to time-varying shifts. The second-order ac 
Stark shift for a two-level system is given by: 



A W P 

ov ac = ~ Re 



1 



1 



A — v + ij A + + 27 J ' 



(5) 



where v is the applied rf frequency, d is the dipole ma- 
trix clement, and A is the energy separation. Because 
the first term is an odd function of detuning, only the 
second so-called Bloch-Siegert term contributes to an ac- 
tual shift of the central frequency of the resonance line- 
shape. For v — A 3> 7 and a transition near satu- 
ration for which d 2 E 2 /j 2 ~ 1, the corresponding shift 



is 8v ac ~ 7 2 /(4A). Hence for 7 = 20 kHz and transi- 
tion frequencies ~ 3 — 1000 MHz, the shift varies from 
~ 0.1 — 30 Hz. For the 3.1 MHz-transition, for which the 
shift is largest, in order to achieve a sensitivity to fre- 
quency shifts of a few mHz, the amplitude stability must 
be better than 10 . However, this requirement becomes 
much less stringent for higher frequency transitions. For 
v ~ 1 GHz, only a modest control at a level of a few 
percent is required. 

We now consider frequency shifts due to all other levels 
on the transition frequency between levels m and n: 



4 V A m? - - v A Tl 



d nk E ( \ 



4 V A nk - v A nk + v 



(6) 



where d m j is the dipole matrix element and A m j = 
E m — Ej is the energy separation between levels m and j. 
The transition widths have been ignored because these 
levels are off-resonant. For even isotopes, this shift is 
~ 0.3 Hz and is mostly determined by levels which are 
> 2000 cm" 1 away with dipole matrix elements ~ 1 
eao. A comparable shift occurs for odd isotopes, assum- 
ing a particular choice of hyperfine transitions for which 
A mj — v~l GHz and d mj — 4 kHz/ (V/cm). Thus, tak- 
ing into account the shifts due to other levels does not 
lead to more stringent requirements on the amplitude 
stability of the rf electric field beyond those obtained in 
the two-level approximation. 



C. Stray magnetic fields 

The residual magnetic field (B) can be controlled to 
~ 1 fiG in the magnetically shielded interaction region. 
This corresponds to Zeeman shifts gp/J-oB ~ 1 Hz, where 
fiQ is the Bohr magneton and gp is the Lande g-factor 
(Table HID . If the Zeeman sublevels are equally popu- 
lated but still unresolved, then this field does not shift 
the central frequency of a resonance lineshape but rather 
only broadens the lineshape with a corresponding width 
of ~ F(gpB — 9FA)fM>B w 1 Hz. However, an imbalance 
of sublevel populations can lead to asymmetric broaden- 
ing of the resonance lineshape, causing an apparent shift 
in the central frequency. In the experimental geometry 
discussed (Fig. |3J|, such an imbalance may be caused by 
residual circular polarization coupled with misalignments 
of the propagation direction of the linearly polarized light 
beams used in the first and second step of the population. 
Because residual circular polarization can be controlled 
to the level of ~ 10~ 5 using standard polarimetric tech- 
niques, the shift of the resonance frequency can be made 
< 10" 5 Hz. 



TABLE III: g-values for the hyperfine levels of states A and 
B ( gjA = 1.21 and g JB = 1.367 pi). 
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FIG. 5: Electric field frequency shifts due to the Millman 
effect. 



E. Black-body radiation 



D. Millman effect in electric resonance 



Misalignments of the atomic beam and the geometry 
of the electrodes can cause the direction of the rf E- 
field to rotate, as seen by a moving atom. This may 
lead to frequency shifts [23j analogous to the Millman 
effect encountered in magnetic resonance |24[- Let Q be 
the frequency of this apparent E-field rotation. Because 
an oscillating field can be decomposed into two counter- 
rotating components, one of these components is shifted 
by +fl and the other by —ft. 

Using the measured homogeneity of our electric field, 
we estimate £1 < 50 Hz. However, for a resonance line- 
shape with unresolved sublevels, the central frequency 
remains largely unaffected. This is not true in the case 
of magnetic resonance (where transitions occur between 
magnetically split sublevels of the same level) and can be 
explained as follows: In the basis in which the quantiza- 
tion axis is perpendicular to the electric field, there are 
only cr + and cr_ transitions. The Millman effect causes 
shifts in the frequencies of these transitions as illustrated 
in Fig.[3]for a A J = transition on resonance (where we 
have assumed that the apparent electric field rotation is 
in the same sense as the <r+ component). If the sublevels 
of the initial state are equally populated and Q <C 7, then 
these shifts only lead to a broadening of the lineshape. 
However, if there is an imbalance of sublevel populations 
the broadening is asymmetric, which can cause an ap- 
parent shift in the central frequency by an amount ~ £f2, 
where £ is the degree of atomic orientation. As men- 
tioned earlier, residual circular polarization in the lin- 
early polarized light, used in the first and second step of 
the population, may induce atomic orientation, but can 
be controlled to the level of ~ 10~ 5 , which is sufficient 
for the desired sensitivity to a. Furthermore, it should 
be noted that it is not the magnitude of the frequency 
shift but rather the stability that is important. 



Black-body radiation (BBR) can cause ac Stark 
shifts ■ The rms value of the black-body radiation 
electric field is 



(E 2 ) = (8.3 V/cm) 2 [T(K)/300 K} 4 



(7) 



A shift in the transition frequency arises from the differ- 
ence in the ac Stark shifts experienced by levels A and B. 
We can give a rough estimate of this Stark shift (5vbb) 
for one of these levels using Eq. ©: 



8v B B 



d 2 (E 2 



4(A - v BB ) ' 



(8) 



where A is a characteristic atomic energy scale, vbb is a 
characteristic frequency of room-temperature BBR, and 
d = 1 eao is a typical optical transition dipole moment. 
Due to cancellations of contributions from nearby energy 
levels above and below the A and B levels, we assume that 
the net shift comes mainly from levels with large energy 
separations: vbb <C A ^ 10 14 Hz, gives an estimate of 
8u BB < 0.1 Hz. 

The interaction region is illuminated by much hot- 
ter BBR from the atomic oven (T~ 1500 K) which is 
^15 cm away. However, the effect is smaller due to the 
decreased solid angle: 



8v 



^ BS [T(K)/300 K] 4 

\ 47T J 
4 / -, 2 

1 cm 



0.1 Hz 



/1500 K 



V 300 K 
= 0.02 Hz, 



4-7T 15 2 cm 2 



(9) 



where we have assumed that an atom in the interaction 
region sees 1 cm 2 of hot surface. 

The BBR shifts are large compared to the desired level 
of sensitivity and are difficult to eliminate. However, 
these shifts can be kept constant to the desired level by 
stabilizing the temperature of the interaction region (to 
~ 2 °C) and the oven (to ~ 30 °C). For the oven, pos- 
sible changes in surface emissivity will be a concern and 
requires further investigation. 
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F. Collisional shifts 

Because the typical cross-sections of the atomic col- 
lisions arc on the order of 1CP 14 cm 2 , most atoms in 
the atomic beam do not experience collisions with other 
atoms before they encounter the back wall of the inter- 
action region. Collisional shifts of ~ 1 — 10 MHz/Torr 
(corresponding to the above-mentioned collisional cross- 
section values) are typically found for atomic transitions 
when the width is measured as a function of the atomic 
pressure. A possible effect of the rare collisions in the 
beam can be estimated by extrapolating these numbers 
to the low pressure values corresponding to the beam 
environment. Assuming a 10 MHz/Torr shift due to col- 
lisions with residual gas, in order to keep time-dependent 
shifts below a few mHz, the residual gas pressure must 
be stable at a level ~ 10~ 10 Torr. This can be achieved 
with standard ultra- high- vacuum equipment. The pres- 
sure stability and composition of the residual gas will be 
monitored with a residual gas analyzer. 

Note that we may find that collisional shifts are, in 
fact, much smaller than the above estimate because the 
A-B transition is nominally between f and d inner-shell 
electrons [2(| . Another important effect that may reduce 
the effect of collisions on the measurement of the vari- 
ation of a is that it is likely that an atom experiencing 
a collision will be quenched to a different atomic state, 
and, thus, will not produce a 564-nm fluorescent photon, 
and will therefore avoid detection. 

A similar consideration applies to collisions between 
Dy atoms themselves. In the absence of collisional 
quenching, the worst-case estimate is that the present 
Dy density (~ 10 10 atoms/cm 3 ) could give a shift of up 
to ~ 2 Hz. This effect will be investigated by measuring 
transition frequencies as a function of the atomic beam 
intensity and, if needed, the Dy beam intensity can be 
monitored and stabilized during measurement. Measur- 
ing a transition frequency at several different Dy-beam 
intensities will allow extrapolation to the collision-free 
value of the frequency. Technically, it is possible to vary 
the density of atomic in the atomic beam without sig- 
nificantly affecting other conditions of the experiment. 
For example, black-body radiation intensity in the inter- 
action region can be maintained constant by separately 
controlling the front and back temperatures of the Dy 
oven. The BBR intensity is determined by the tempera- 
ture of the front part of the oven near the nozzle, while 
the flux of atoms is determined by the lower temperature 
towards the back of the oven. 

Another concern is a possible collisional wall-shift re- 
sulting from a fraction of the atoms reentering the in- 
teraction region upon reflection from its back wall. This 
effect is suppressed by a number of small factors: reflec- 
tion probability, solid angle for entering the interaction 
region, the probability for an atom to avoid quenching 
in the wall collision, etc. If necessary, cryo-cooling of the 
interaction-region walls may be employed to reduce the 
probability for an atom to bounce from it. 




Ml 



A B 

FIG. 6: Simultaneous measurement of an Ml transition fre- 
quency for a nonzero nuclear spin isotope facilitates the con- 
trol of systematics. 



G. Doppler Shift 

The effect of the first-order Doppler shift is estimated 
to be small. To see this, we model the electric field 
(Fig. as a standing wave constructed from two travel- 
ing waves counter-propagating co- linearly with the direc- 
tion of the atomic beam. Due to ohmic losses inside the 
conductors, the amplitude of each wave gets attenuated. 
A difference in the intensity of these waves at the location 
of the atoms leads to asymmetric Doppler broadening of 
the lineshape, and thus an apparent shift. We estimate 
this intensity difference by first considering the B-field 
induced by the time-varying E-field. This B-field pene- 
trates inside the metal to within a skin-depth, inducing 
currents from which the power loss can be readily cal- 
culated [27|. This simple model gives a fractional power 
difference oc z/ 5 / 2 . Since the Doppler shift for a given 
traveling- wave power is oc ^, the shift is oc v 1 1 2 . For a 
transition frequency of 1 GHz, we estimate the asymme- 
try to be < 10~ 4 and the corresponding shift is ~ 0.2 Hz. 
One can imagine a factor of ~ 10 suppression if rf power 
is fed to the plates in a symmetric fashion. Thus, it is 
only required that this shift be stable to ~ 10% in order 
to achieve a sensitivity of a few mHz. 

In addition to the first-order Doppler shift, we also 
consider the second-order effect. Depending upon the rf 
transition, the second order Doppler shift is ~ 10 -3 — 
10 -6 Hz, which is sufficiently small. 



H. Techniques to control systematics 

A powerful method to detect and eliminate possible 
sources of systematic shifts common to both levels is to 
simultaneously measure the transition frequency between 
hyperfine levels of a given parity state. The reason is 
that, because the levels involved have the same relativis- 
tic corrections, this frequency is insensitive to variations 
of a. One possible scheme is to excite an Ml transition, 
e. g., as shown in Fig. whose frequency can be mon- 
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FIG. 7: Simultaneous measurement of two (or more) El tran- 
sition frequencies to facilitate the control of systematics. This 
method works with the two zero-spin isotopes, 162 Dy and 
164 Dy. 

itored by looking for disappearance in the fluorescence 
from level A for a fixed El transition frequency. Alter- 
natively, we can utilize another El transition as shown 
in Fig. In this scenario, the effect a variation is twice 
as large in the sum of the two frequencies, while the dif- 
ference is insensitive to a variations. 

Furthermore, one can compare El transitions for the 
two abundant isotopes with zero nuclear spin ( 162 Dy and 



[1] J. Uzan, Rev. Mod. Phys. 75, 403 (2003). 
[2] W.J. Marciano, Phys. Rev. Lett. 52, 489 (1984). 
[3] J.D. Barrow, Phys. Rev. D 35, 1805 (1987). 
[4] T. Damour and A.M. Polyakov, Nucl. Phys. B423, 532 
(1994). 

[5] T. Damour, F. Piazza, and G. Veneziano, Phys. Rev. 

Lett. 89, 081601 (2002). 
[6] J. K. Webb, et al. , Phys. Rev. Lett. 87, 091301 (2001). 
[7] A.I. Shlyakhter, Nature 264, 340 (1976). 
[8] T. Damour and F. Dyson, Nucl. Phys. B480, 37 (1996). 
[9] Y. Fujii, et al. , Nucl. Phys. B573, 377 (2000). 
[10] S. K. Lamoreaux, (2003) e-print physics/0309048 at 

http: / /arxiv.org/ 
[11] J. D. Prestage, R. L. Tjoelker, and L. Maleki, Phys. Rev. 

Lett. 74, 3511 (1995). 
[12] H. Marion, et al. , Phys. Rev. Lett. 90, 150801 (2003). 
[13] S. Bize, et al., Phys. Rev. Lett. 90, 150802 (2003). 
[14] V. A. Dzuba, V. V. Flambaum, and J. K. Webb, Phys. 

Rev. Lett. 82, 888 (1999). 
[15] V. A. Dzuba, V. V. Flambaum, and J. K. Webb, Phys. 

Rev. A 59, 230 (1999). 
[16] V. A. Dzuba, V. V. Flambaum, and M. V. Marchenko, 



164 Dy). The counting rate is significantly higher and the 
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V. DISCUSSION 

In summary, rf El transitions in Dy provide an attrac- 
tive system in which to test the temporal variation of 
a. The frequencies of these transitions can be directly 
counted. For a limit of a/ a < 10 _15 /yr, the shift was 
calculated to be ~ 2 Hz/yr. At present, a statistical 
sensitivity of 0.6 Hz in one second of integration time is 
achievable. Knowledge of systematic effects is critical to 
this experiment. Preliminary analysis shows that it may 
be possible to control them at a level corresponding to 
\a/a\ ~ 10~ 18 /yr, a level of sensitivity that would rival 
that of the most stringent observational limit set by the 
Oklo natural reactor. 
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